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ABSTRACT 

We present preliminary results of a new measurement of the inclusive b quark fragmen- 
tation function in Z° decays using a novel kinematic B hadron energy reconstruction 
technique. The measurement is performed using 150,000 hadronic events recorded 
in the SLD experiment at SLAC between 1996 and 1997. The small and stable SLC 
beam spot and the CCD-based vertex detector are used to reconstruct topological B- 
decay vertices with high efficiency and purity, and to provide precise measurements of 
the kinematic quantities used in this technique. We measure the B energy with good 
efficiency and resolution over the full kinematic range. We compare the measured scaled 
B hadron energy distribution with several functional forms of the B hadron energy dis- 
tribution and predictions of several models of b quark fragmentation. Several functions 
are excluded by the data. The average scaled energy of the weakly decaying B hadron is 
measured to be = 0.714 ± 0.005 (stat) ± 0.007 (syst) ± 0.002 (model) (preliminary). 
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1. Introduction 



The production of heavy hadrons (H) in e^e~ annihilation provides a laboratory for 
the study of heavy-quark (Q) jet fragmentation. This is commonly characterised in 
terms of the observable xh = 'iEu / where Eh is the energy of a B or D hadron 
containing a 6 or c quark, respectively, and ^/s is the cm. energy. In contrast to 
light-quark jet fragmentation one expects [0 the distribution of xh, D{xh), to peak at 
an xj^- value significantly above 0. Since the hadronisation process is intrinsically non- 
perturbative D{xh) cannot be calculated directly using perturbative Quantum Chromo- 
dynamics (QCD). However, the distribution of the closely- related variable xq = 2EqI ^ 
can be calculated perturbatively 0, |, Q and related, via model-dependent assumptions, 
to the observable quantity D{xh)', a number of such models of heavy quark fragmenta- 
tion have been proposed [|, ||, 0]. Measurements of D{xh) thus serve to constrain both 
perturbative QCD and the model predictions. Furthermore, the measurement of D{xh) 
at different cm. energies can be used to test QCD evolution, and comparison of D{xb) 
with D{x£)) can be used to test heavy quark symmetry P, |^. Finally, the uncertainty on 
the forms of D^xn) and D{xb) must be taken into account in studies of the production 



and decay of heavy quarks, see eg. [T^ ; more accurate measurements of these forms will 
allow increased precision in tests of the electroweak heavy-quark sector. 

We consider the measurement of the B hadron scaled energy distribution D{xb) 



in Z° decays. Earlier studies 0] used the momentum spectrum of the lepton from 



semi-leptonic B decays to constrain the mean value < xb > and found it to be approx- 
imately 0.70; this is in agreement with the results of similar studies at ^/s = 29 and 35 
GeV |12|. In more recent analyses [13, 14, 15] the scaled energy distribution D[xb) has 



been measured by reconstructing B hadrons via their B D/X decay mode. In this 
case the reconstruction efficiency is intrinsically low due to the small branching ratio 
for B hadrons to decay into the high-momentum leptons used in the tag. Also, the 
reconstruction of the B hadron energy using calorimeter information usually has poor 
resolution for low B energy, resulting in poor sensitivity to the shape of the distribution 
at low energy. 

Here we describe the preliminary results of a new method for reconstructing B hadron 
decays, and the B energy, inclusively, using only charged tracks, in the SLD experiment 
at SLAC. We use the upgraded CCD vertex detector, installed in 1996, to reconstruct 
5-decay vertices with high efficiency and purity. Combined with the micron-size SLC 
interaction point (IP), precise vertexing allows us to reconstruct accurately the B flight 
direction and hence the transverse momentum of tracks associated with the vertex with 
respect to this direction. Using the transverse momentum and the total invariant mass 
of the associated tracks, an upper limit on the mass of the missing particles is found for 
each reconstructed 5-decay vertex, and is used to solve for the longitudinal momentum 
of the missing particles, and hence for the energy of the B hadron. In order to improve 
the B sample purity and the reconstructed B hadron energy resolution, B vertices with 
low missing mass are selected. The method is described in Section 3. In Section 4 
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we compare the B energy distribution with predictions of heavy quark fragmentation 
models. We also test several functional forms of B hadron energy distributions. In 
Section 5, we unfolded the B hadron energy distribution. In Section 6, we discuss the 
systematic errors. In Section 7 we summarize the results. 



2. Apparatus and Hadronic Event Selection 

This analysis is based on roughly 150,000 hadronic events produced in e'^e~ annihilations 
at a mean center-of-mass energy of ^/s = 91.28 GeV at the SLAG Linear Collider 
(SLC), and recorded in the SLC Large Detector (SLD) in 1996 and 1997. A general 
description of the SLD can be found elsewhere |]16|. The trigger and initial selection 
criteria for hadronic decays are described in Ref. |]T^. This analysis used charged 
tracks measured in the Central Drift Chamber (CDC) and in the upgraded Vertex 
Detector (VXDS) |jT9|. Momentum measurement is provided by a uniform axial magnetic 
field of 0.6T. The CDC and VXDS give a momentum resolution of <Jpj^/p± = 0.01 © 
0.0026j9_L, where p± is the track momentum transverse to the beam axis in GeV/c. In the 
plane normal to the beamline the centroid of the micron-sized SLC IP is reconstructed 
from tracks in sets of approximately thirty sequential hadronic decays to a precision 
of a'^'^ ^ 7±2 /im (1996) and a^^ ^ 4±2 fim (1997). The IP position along the beam axis 
is determined event by event using charged tracks with a resolution of ~ S5 /xm (1996) 
and ~ SO fim (1997). Including the uncertainty on the IP position, the resolution 
on the charged-track impact parameter [d) projected in the plane perpendicular to the 
beamline is a^"^ = 14©SS/(psin3/2 e) /xm (1996) and a^/ = ll©33/(psin3/2 ^) (1997), 
and the resolution in the plane containing the beam axis is cr^ = 27©SS/(psin^^^ ^) /im 
(1996) and crj = 24©SS/(psin^'^^ 6*) /im (1997), where 6 is the track polar angle with 
respect to the beamline. The event thrust axis |]20| is calculated using energy clusters 



measured in the Liquid Argon Calorimeter [21 



A set of cuts is applied to the data to select well-measured tracks and events well 
contained within the detector acceptance. Charged tracks are required to have a distance 
of closest approach transverse to the beam axis within 5 cm, and within 10 cm along 
the axis from the measured IP, as well as | cos 6^1 < 0.80, and p± > 0.15 GeV/c. Events 
are required to have a minimum of seven such tracks, a thrust axis polar angle w.r.t. 
the beamline, 9^, within | cos 6*^1 < 0.71, and a charged visible energy E^s of at least 
20 GeV, which is calculated from the selected tracks assigned the charged pion mass. 
The efficiency for selecting a well-contained Z^ QQio) event is estimated to be above 
96% independent of quark flavor. The selected sample comprised 111,569 events, with 
an estimated 0.10 ± 0.05% background contribution dominated by Z° t^t~ events. 

For the purpose of estimating the efficiency and purity of the B hadron selection 
procedure we made use of a detailed Monte Carlo (MC) simulation of the detector. 



The JETSET 7.4 p2[ event generator is used, with parameter values tuned to hadronic 



e^e annihilation data ||2^, combined with a simulation of B hadron decays tuned ||24 



to T(4S') data and a simulation of the SLD based on GEANT S.21 Inclusive 
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distributions of single-particle and event-topology observables in hadronic events are 
found to be well described by the simulation [1^. Uncertainties in the simulation are 
taken into account in the systematic errors (Section ^). 



3. B Hadron Selection and Energy Measurement 
3.1. B Hadron Selection 

The B sample for this analysis is selected using a topological vertexing technique based 
on the detection and measurement of charged tracks, which is described in detail in 
Ref. 1 2^. Each hadronic event is divided into two hemispheres by a plane perpendicular 



to the thrust axis. In each hemisphere the topological vertexing algorithm is applied to 
the set of ' quality' tracks having (i) at least 23 hits in the CDC and 2 hits in VXDS; 
(ii) a combined CDC and VXDS track fit quality of x^/^dof <8; (iii) a momentum in 
the range 0.25< p <55 GeV/c, (iv) an impact parameter of less than O.S cm in the r(j) 
plane, and less than 1.5 cm along the z axis; (v) a transverse impact parameter error no 
larger than 250 /xm. 

Vertices consistent with photon conversions or and A° decays are discarded. In 
hemispheres containing at least one found vertex the vertex furthest from the IP is 
retained as the 'seed' vertex. Those events are retained which contain a seed vertex 
separated from the IP by between 0.1 cm and 2.S cm. The lower bound reduces con- 
tamination from non-i?-decay tracks and backgrounds from light-flavor events, and the 
upper bound reduces the background from particle interactions with the beam pipe. 

For each hemisphere containing an accepted seed vertex, a vertex axis is formed by 
the straight line joining the IP to the seed vertex, which is located at a distance D from 
the IP. For each quality track not directly associated with the vertex, the distance of 
closest approach to the vertex axis, T, and the distance from the IP along the vertex 
axis to the point of closest approach, L, are calculated. Tracks satisfying T< 1 mm and 
L/D> O.S are added to the vertex. These T and L cuts are chosen to minimize false 
track associations to the seed vertex, since typically the addition of a false track has a 
much greater kinematic effect than the omission of a genuine S-decay track, and hence 
has more effect on the reconstructed B hadron energy resolution. Our Monte Carlo 
studies show that, on average, this procedure attaches 0.85 tracks to each seed vertex, 
91.9% of the tracks from tagged true B decays are associated with the resulting vertices, 
and 98.0% of the vertex tracks are from true B decays. 

The large masses of the B hadrons relative to light-flavor hadrons make it possible to 
distinguish B hadron decay vertices from those vertices found in events of light flavors 
using the vertex invariant mass, M. However, due to the missing particles, which are 
mainly neutrals, M cannot be fully determined. In the rest frame of the decaying hadron, 
M can be written as 



M = JM^, + + PI, + J Mi + P2 + p2 (1) 
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where Mch and Mq are the total invariant masses of the set of vertex-associated tracks 
and the set of missing particles, respectively. Pt is the total charged track momentum 
transverse to the B flight direction, which is identical to the transverse momentum of 
the set of missing particles by momentum conservation. Pchi and Pqi are the respective 
momenta along the B flight direction. In the B rest frame, Pcm = Pqi- Using the set 
of vertex- associated charged tracks, we calculate the total momentum vector Pch-, the 
total energy E^h and the invariant mass M^.^, assuming the charged pion mass for each 
track. The B hadron flight direction (the line joining the IP and the B vertex. The 
lower bound for the mass of the decaying hadron, the 'Pf-corrected vertex mass'. 



Mp, = ^Ml^ + P? + \Pt\ (2) 

is used as the variable for selecting B hadrons. The majority of non-i? vertices have Mpt 
less than 2.0 GeV/c^. However, occasionally the measured Pt may fluctuate to a much 
larger value than the true P^, causing some charm vertices to have a Mpt larger than 2.0 
GeV/c^. To reduce this contamination, we calculate the 'minimum Pt by allowing the 
locations of the IP and the vertex to float to any pair of locations within the respective 
one sigma error-ellipsoids. We substitute the minimum Pt in Equation and use the 
modified Mpt as our variable for selecting B hadrons [ P7[] . 

Figure |I| shows the distribution of the Mpt for the 32,492 hemispheres in the data 
sample with a found secondary vertex, and the corresponding simulated distribution 
(histogram). B hadron candidates are selected by requiring Mpt > 2.0 GeV/c^. We 
further required Mpt < 2 x M^h to reduce the contamination from fake vertices in 



light quark events [^. A total of 19,604 hemispheres are selected, with an estimated 
efficiency for selecting a true P-hemisphere of 40.1%, and a sample purity of 98.2%. 
The contributions from light-flavor events in the sample are 0.15% for primary u,d and 
s events and 1.6% for c events. 



3.2. B Hadron Energy Measurement 

The energy of each B hadron, Ep, can be expressed as the sum of the reconstructed- 
vertex energy, Ech, and the energy of those particles not associated with the vertex, 

Eq. 

We can write Eq as 

Po' = + p2 + p2 (3) 

The two unknowns, Mq and Pqi, must be found in order to obtain Pq. One kinematic 
constraint can be obtained by imposing the P hadron mass on the vertex, Mp = Ep—Pp, 
where Pp = PcM + Pqi is the total momentum of the P hadron, and Pchi is the momentum 
component of the vertex-associated tracks along the vertex axis. From Equation (|I]) we 
derive the following inequality, 

^jM^, + Pt' + ^Mi + Pt'<Mp, (4) 
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where equality holds in the limit where both Pqi and Pchi vanish in the B hadron rest 
frame. Equation (|^) effectively sets an upper bound on Mo, and a lower bound is given 
by zero: 

< < Mo^^,, (5) 

where 



MLa. = M's- 2Mb^MI + P2 + M^,. (6) 

Since Mq is bounded from both above and below, we expect to obtain a good estimate 
of Mq, and therefore of the B hadron energy, when Mq^^^ is small. 

We have used our simulation to study this issue. Assuming Mb = 5.28 GeV/c^, 
the true value of Mq tends to cluster near its maximum value Momax- Figure |^ shows 
the relative deviation of Momax from Mgtrue for all B hadrons. Although approximately 
20% of the B hadrons are 5° and A;, which have larger masses, the values of M^^ax 
obtained using Mb=5.28 GeV/c^ in Equation (|^) are typically within about 10% of 
Mq. The distribution of the reconstructed Mq^„^ for vertices in the selected B hadron 
sample is shown in Figure ^. The simulation indicates that the non-66 background is 
concentrated at high Mq^^^; this because most of the light flavor vertices have small 
Mpt and therefore, due to the strong negative correlation between Mpt and Mo^ax, 
large Mq^ux- The negative tail in Figure |^ is an effect of detector resolution, and the 
Monte Carlo simulation shows good agreement with the data. 

Because Mq peaks near Mo^nax, we set Mq = Mq^^^ if Mq^^^ >0, and Mq = if 
^Omax <0- We then calculate Pqi: 



p _ Ml-{Ml, + Pf)-{Ml + Ph ^ 
Poi - 2(MlTif^ 



and hence Eq (Equation (^)). We then divide the reconstructed B hadron energy, 
£■^=0 ^ _|_ j^^^^ ^YiQ beam energy, Ebeam = to obtain the reconstructed scaled 

B hadron energy, x^'p" = E'^^'/Ebeam- 

The resolution of x^*^ depends on both Mq^„^ and the true xb, x^^'^. Vertices in 
the negative tail of the Mq^„^ distribution that have Mq^^^. < —1.0{GeV/c^)'^ are often 
poorly reconstructed and are not used in further analysis. Vertices with small values of 
l^omaxl typically reconstructed with better resolution and an upper cut on Mq^^^ 
is hence applied. For an x^-independent cut, the efficiency for selecting B hadrons is 
roughly linear in x^"^. In order to obtain an approximately x^-independent selection 
efficiency we choose the following upper cut: 

Mq2^,, < {1.1 + 0.006(E,,,^ - ES^-^) + S-bexpl-iE^B^" ' 5.5)/3.5]}' , (8) 

where the two terms that depend on the reconstructed energy E^'^ increase the efficiency 
at lower B hadron energy. 

Only about 0.7% of the selected vertices are from light-flavor events, but they are 
concentrated in the lowest energy bin. To further remove this background, a vertex is 



6 



required to contain at least 3 quality tracks with a normalized impact parameter greater 
than 2. This eliminates almost all of the uds-event background and cuts the charm 
background by about 20% overall and 43% in the few lowest energy bins. This cut 
helps to reduce the dependence of the reconstructed B hadron energy distribution on 
the light flavor simulation in the low energy region, which is a key step towards finding 
the correct shape of the B hadron energy distribution at low energies. Figure | shows 
the distribution of M^^^^ after all these cuts; the data and Monte Carlo simulation are 
in good agreement. 

A total of 1920 vertices in the data for 1996-97 satisfy all these selection cuts. The 
overall efficiency for selecting B hadrons is 3.9% and the estimated B hadron purity is 
99.5%. The efficiency as a function of x^"*^ is shown in Figure ^. The dependence is 
rather weak except for the lowest xb region; the efficiency is substantial, about 1.7% 
even just above the kinematic threshold for B energy. 

We examine the i?-energy resolution of this technique. The distribution of the 
normalized difference between the true and reconstructed B hadron energies, (x^'^ — 
for Monte Carlo events, is fitted by a double Gaussian, resulting in a core 
width (the width of the narrower Gaussian) of 10.4% and a tail width (the width of the 
wider Gaussian) of 23.6% with a core fraction of 83%. Figure ^ shows the core and tail 
widths as a function of x%'^'^. In order to compare the widths from different xb bins, 
we fix the ratio between core and tail fractions to that obtained in the overall fit above. 
The x^-dependence of the resolution is weak, indicating that the absolute resolution 
on Xb, x'^b'^ — x%'^'^, is very good at low B energy, which is an advantage of this energy 
reconstruction technique. 

Figure |^ shows the distribution of the reconstructed scaled B hadron energy for 
the data, D'^"*"(a;^'^), and for the Monte Carlo simulation, D^'^(a;^'^). The small non- 
66 background, the high B selection efficiency over the full kinematic coverage, and the 
good energy resolution combine to give a much improved sensitivity of the data to the 
underlying true shape of the B energy distribution (see next section). 

The event generator used in our simulation is based on a perturbative QCD 'par- 
ton shower' for production of quarks and gluons, together with the phenomenological 
Peterson function |^ (Table Q) to account for the fragmentation of b and c quarks into 
B and D hadrons, respectively, within the iterative Lund string hadronisation mecha- 
nism p^; this simulation yields a 'generator-level' primary B hadron energy distribution 



with < Xb > = 0.6930. It is apparent that this simulation does not reproduce the data 
well (Figure |^); the for the comparison is 62 for 16 binsQ 

The distribution of the non-66 background, S'(a;^'^), is also shown in Figure ^ The 
background is subtracted bin- by-bin from the Z}'^"*'*(xg^) before we proceed to test var- 
ious fragmentation models. 



*We used a value of the Peterson function parameter et = 0.006 pg 
^We exclude several bins with very few events in the comparison, 
details. 



For details see Section 4.1 for 
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4. The Shape of the B Hadron Energy Distribution 



Given the raw reconstructed B energy distribution in the data shown in Figure 0, there 
are several ways of estimating the true underlying B energy distribution. Here we take 
two approaches, each described in a subsection. 

In the first part, we test several b fragmentation models, f{z,f5) embedded within 
Monte Carlo generators, where z is an internal, experimentally inaccessible variable, 
corresponding roughly to the fraction of the momentum of the fragmentating b quark 
carried by the resulting B hadron, and f3 is the set of parameters associated with the 
model in question. In the second part, we test several functional forms for the distribu- 
tion of xb itself, /(xb. A), where A represents the set of parameters associated with each 
functional form. 



4.1. Tests of b Quark Fragmentation Models f{z^(3) 

We first consider testing models of b quark fragmentation. Since the fragmentation func- 
tions for various models are usually functions of an experimentally inaccessible variable 
(e.g. z = {E+p\\)h/ {E+p\\)q or z = VWhIvWq )i it is necessary to use a Monte Carlo gen- 
erator to generate events according to a given input heavy quark fragmentation function 
f{z,(3), where P represents the set of parameters. 

We consider the phenonmenological models of the Lund group 0, Bowler 0, Pe- 
terson et al. 1^ and Kartvelishvili et al. ||2^. We also consider the perturbative QCD 



calculations of Braaten et a/.(BCFY) [|], and of Collins and Spiller (CS) [Q. We use 
the JETSET parton shower Monte Carlo and each fragmentation model in question 
to generate the simulated events without detector simulation. Table |l| contains a list of 
the models. In addition, we test the UCLA fragmentation model which with fixed 



parameters. For b fragmentation, we also test the HERWIG ||3T]] event-generator using 
both possible values of the parameter cldir = and 1. 

In order to make a consistent comparison of each model with the data we adopt the 
following procedure. For each model, starting values of the arbitrary parameters, (3, are 
assigned and the corresponding fragmentation function f{z,P) is used along with the 
JETSET Monte Carlo to produce the corresponding scaled primary B hadron energy 
distribution, D^''^{x^^^'^) in the MC-generated bb event sample, before simulation of the 
detector. Then each simulated B hadron is weighted according to its true B hadron 
energy, x^**^; the weight is determined by the ratio of the generated B hadron energy 
distribution, D^'^^ {x%^^), to that of our default simulation jjd'^f^-^^ti^^truey After simu- 
lation of the detector, application of the analysis cuts and background subtraction, the 
resulting weighted distribution of reconstructed B hadron energies, D^^'-^{x'''^'^), is then 
compared with the background-subtracted data distribution and the value, defined 
as 

X= . ^ PL-^ (9) 
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Model 




Reference 


BCFY 
Bowler 

CS 
Kart. 
Lund 
Peterson 


n ^^].~^K i(i[3 + E- i( zYfAr)] 
[1 - (1 - r)z]^^ ^'-^^ ' -"^ 

^(i+r!bml)(l z)'^ea:p( hml/ z) 
In 1 ^-2 


i 

29 

i 
i 



Table 1: h quark fragmentation models used in comparison with the data. For the BCFY 
model, /i(r) = 3(3 - 4r), = 12 - 2?>r + 2Qr'^ , f^{r) = (1 - r)(9 - llr + 12r2), 
and/4(r) = 3(1 - r)2(l - r + r^). 



is calculated, where is the number of bins to be used in the comparison, 
number of entries in bin i in the data distribution, and N^'-^ is the number of entries 
in bin i in the simulated distributionl^. cTj is the statistical error on the deviation of the 
observed number of entries for the data from the expected number of entries in bin i, 
which can be expressed as 

= (r^y , (10) 



where y^J tN^'^^ j is the expected statistical variance on the observed data number of 

entries in bin z, assuming the model being tested is correct, and (^r^J^^f^^ is the 

statistical variance on the expected number of entries in bin i. Since the x^-test is not 
a statistically effective test for bins with a very small number of entries, the third, the 
fourth, and the last three bins in Figure |^ are excluded from the comparison. 

We vary the values of the set of parameters j3 and repeat the above procedure. The 
minimum is found by scanning through the input parameter space, yielding a set 
of parameters which give an optimal description of the reconstructed data by the frag- 
mentation model in question. Each of the nine plots in Figure shows the background- 
subtracted distribution of reconstructed B hadron energy for the data (points) and the 
respective B energy distribution (histogram) resulting either from the optimised input 
fragmentation function f{z) embedded within the JETSET parton shower simulation, 
or from the predictions of the HERWIG event-generator and the UCLA fragmentation 
model. Data points excluded from the fit are represented in Figure | by open circles. 
Table |^ lists the results of the comparisons. 

is the factor by which the total number of entries in the simulated distribution is scaled to the 
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Model 


xVdof 


Parameters 


{xb) 


JETSET + BCFY 


83/16 


r = 0.085 


0.694 


JETSET + Bowler* 


17/15 


a= 1.5,6 = 1.5, (rfe = 1) 


0.714 


JETSET + Collins and Spiller 


103/16 


efc = 0.003 


0.691 


JETSET + Kartvelishvili* et al 


34/16 


ttfe = 10.4 


0.711 


JETSET + Lund* 


17/15 


a = 2.0,6 = 0.5 


0.712 


JETSET + Peterson et al. 


62/16 


eb = 0.006 


0.697 


HERWIG cldir=0 


460/17 




0.632 


HERWIG cldir=l 


94/17 




0.676 


UCLA* 


25/17 




0.718 



Table 2: Results of fragmentation model tests for JETSET + fragmentation models, the 
HERWIG model and the UCLA model. Minimum x^, number of degrees of freedom, 
coresponding parameter values, and the mean value of the corresponding B energy 
distribution are listed. * indicates used to correct the data in Section |^. 



We conclude that with our resolution and our current data sample, we are able to 
distinguish between several fragmentation models. Within the context of the JETSET 
Monte Carlo, the Lund and Bowler models are consistent with the data with proba- 
bility of 32% for each, the Kartvelishvili model is marginally consistent with the data, 
while the Peterson, the BCFY and the CS models are found to be inconsistent with the 
data. The UCLA model is consistent with the data to a level of 10% probability. 
The HERWIG model with cldir = is confirmed to be much too soft. Using cldir = 1 
results in a substantial improvement, but it is still inconsistent with the data. 



4.2. Tests of Functional Forms /(a:5,A) 

We then consider the more general question of what functional forms of the B energy 
distribution, f{xB, A), can be used as estimates of the true underlying B energy distri- 
bution. In particular, we would like to test a wide variety of functional forms and ask 
how many different forms are consistent with the data. Each consistent functional form 
will add to the list of our estimates of the true underlying B energy distribution. 

For convenience we consider the functional forms of the BCFY, Collins and Spiller, 
Kartvelishvili, Lund, and Peterson groups in the variable xb- In addition we consider ad 
/ioc generalisations of the Peterson function (F), an 8th-order polynomial and a 'power' 
function. These functions are listed in Table Each function vanishes a.t xb = and 
Xb = 1- 

For each functional form, a testing procedure similar to that described in subsec- 
tion [4.1| is applied. The optimised fitting parameters A and the minimum values are 



number of entries in the data distribution; r ~ 1/12. 
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Function 


fixB,X) 


Reference 


F 

8th-order Polynomial 
Power 


{1 + b{l - xb)) c d \-2 

Xb ^ Xb 1-Xb{ 
Xb{1 - Xb){xb - X%){1 + Y!i=iPiXB) 

x%{1-xbY 


Hi 

(see text) 
(see text) 



Table 3: B energy functional forms used in comparison with the data. A polynomial function 
and a power function are included (see text for discussion), is the low kinematic threshold 
for B energy. For BCFY, CS, Kartvelishvili, Lund, Peterson functional forms, see Table 

listed in Table The corresponding D*^'"(x^'^) are compared with the data in Figure]^. 

Two sets of optimised parameters are found for the generalised Peterson function F 
to describe the data. 'Fl', obtained by setting the parameter h (shown in Table ^) to 
infinity, behaves like xb as — >• and (1 — xbY as — >• 1 and yields the best 
probability of 53%; 'F2', obtained by setting b to zero, has a probability of 13%. A 
constrained polynomial of at least 8th-order is needed to obtain a x^ probability greater 
than 0.1%. The Peterson functional form marginally reproduces the data with a x^ 
probability of about 3%. The remaining functional forms are found to be inconsistent 
with our data. The widths of the BCFY and CS functions are too large to describe 
the data; Kartvelishvili, Lund and the 'power' functional form vanish too fast as xb 
approaches zero. We conclude that, within our resolution and with our current data 
sample, we are able to distinguish between some of these functional forms. But most 
importantly, consistent functional forms will help us evaluate the uncertainty on the true 
B energy distribution. 



5. Correction of the B Energy Distribution 

In order to compare our results with those from other experiments and potential future 
theoretical predictions it is necessary to correct the reconstructed scaled B hadron energy 
distribution D'^"*"(x^^) for the effects of non-i? backgrounds, detector acceptance, event 
selection and analysis bias, and initial-state radiation, as well as for bin-to-bin migration 
effects caused by the finite resolution of the detector and the analysis technique. Due to 
the known rapid variation of the yet-unknown true B energy distribution at large xb, 
any correction procedure will necessarily be more or less model-dependent. We choose a 
method that explicitly evaluates this model-dependence and gives a very good estimate 
of the true energy distribution using all of the above models or functional forms that 
are at least marginally consistent with the data. 

We apply a 25 x 25 matrix unfolding procedure to D'^"*"(a;^'^) to obtain an estimate 
of the true distribution L)'^"*"(x^"^): 

jjdata^^true-^ = e-^(a;^™") ■ ^(x*^"^ x'^^) ■ {D^""'^ {x'^") - Six'^")) (11) 
where 5 is a vector representing the background contribution, E is a, matrix to correct 
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Function 


xVdof 


Parameters 


{xb) 


Fl* 


14/15 


c = 


0.838 ±0.018 


0.714±0.005 






d = 


0.022 ±0.002 




F2* 


21/15 


c = 


0.896 ±0.033 


0.717±0.005 






d = 


0.040 ± 0.003 




BCFY 


62/16 


r = 


0.240 ± 0.009 


0.709±0.005 


Collins and Spiller 


75/16 


= 


0.043 ±0.005 


0.711±0.005 


Kartvelishvili et al. 


68/16 




= 4.16 ±0.11 


0.721±0.004 


Lund 


115/15 


a - 


= 2.30 ±0.12 


0.721±0.005 






hm\ 


= 0.50 ±0.07 




Peterson et al* 


28/16 


efe = 


0.036 ±0.002 


0.713±0.005 


Polynomial* 


15/12 


Pi = 


-10.76 ±0.16 


0.709±0.005 






P2 = 


= 45.74 ±0.28 








P3 = 


-93.60 ±0.34 








Pi = 


= 92.01 ±0.37 








P5 = 


-34.53 ±0.27 




Power 


68/15 


a - 


= 4.27 ±0.25 


0.720±0.005 






P- 


= 1.05 ±0.10 





Table 4: Results of the fit of fragmentation functions to the reconstructed B hadron 
energy distribution after background subtraction. The minimum value, the number 
of degrees of freedom, the coresponding parameter values, and the mean value of the 
corresponding B energy distribution are listed. Errors are statistical only. * indicates 
used to correct the data in Section |. 



for bin-to-bin migrations, and e is a vector representing the efficiency for selecting true 
B hadron decays for the analysis. The matrices S", E and e are calculated from our MC 
simulation; the matrix E incorporates a convolution of the input fragmentation function 
with the resolution of the detector. E{i,j) is the number of vertices with x^"*^ in bin i 
and x^^"^ in bin j, normalized by the total number of vertices with x^^"^ in bin j. 

We evaluate the matrix E using the Monte Carlo simulation weighted according to 
an input generator-level true B energy distribution found to be consistent with the data 
in Section |. We have seen that several B energy distributions can reproduce the data. 
We consider in turn each of these eight consistent distributions, using the optimised 
parameters listed in Table |^ and |[ The matrix E is then evaluated by examining the 
population migrations of true B hadrons between bins of the input scaled B energy, 
x^""^, and the reconstructed scaled B energy, x^^'^. Using each D'^'^' {x%'^'^), the data 
distribution D'^^°'{x'^^'^) is then unfolded according to Equation ( 
which is shown for each input fragmentation function in Figure |10|. 

It can be seen that the shapes of ]jdat'^(^x%^'^) differ systematically among the input 
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b quark fragmentation models and the assumed B energy functional forms. These dif- 
ferences are used to assign systematic errors. Figure |Tl| shows the final corrected xb 
distribution D{xb), which is the bin-by-bin average of the eight unfolded distributions, 
where the inner error bar represents the statistical error and the outer error bar is the 
sum in quadrature of the r.m.s. of the eight unfolded distributions and the statistical 
error within each bin. Since two of the eight functions (the Kartvelishvili model and 
the Peterson functional form) are only in marginal agreement with the data, and the 
8th-order polynomial has a slightly unphysical behavior near a;^ = 1, this r.m.s. may be 
considered to be a rather reasonable envelope within which the true xb distribution is 
most likely to vary. The model dependence for this analysis is significantly smaller than 
those of previous direct B energy measurements, indicating the enhanced sensitivity of 
our data to the underlying true energy distribution. 



6. Systematic Errors 

We have considered sources of systematic uncertainty that potentially affect our mea- 
surement of the B hadron energy distribution. These may be divided into uncertainties 
in modelling the detector and uncertainties on experimental measurements serving as 
input parameters to the underlying physics modelling. For these studies our standard 
simulation, employing the Peterson fragmentation function, is used. 

For each source of systematic error, the Monte Carlo distribution D^'-' {x^^'^'^) is re- 
weighted and then the resulting Monte Carlo reconstructed distribution D'^'"{x^^'^) is 
compared with the data by repeating the fitting and unfolding procedures 

described in Section 4 and 5. The differences in both the shape and the mean value 
of the x^"*^ distribution relative to the standard procedure with nominal values of pa- 
rameters are considered. Due to the strong dependence of our energy reconstruction 
technique on charged tracks, the dominant systematic error is due to the discrepancy 
in the charged track transverse momentum resolution between the Monte Carlo and the 
data. We evaluate this conservatively by taking the full difference between the nominal 
results and results using a resolution-corrected Monte Carlo event sample. The differ- 
ence between the measured and simulated charged track multiplicity as a function of 
cos^ and momentum is attributed to an unsimulated tracking inefficiency correction. 
We use a random track-tossing procedure to evaluate the difference in our results. 

A large number of measured quantities relating to the production and decay of charm 
and bottom hadrons are used as input to our simulation. In bb events we have considered 
the uncertainties on: the branching fraction for — >■ bb; the rates of production of B^, 
B^ and B^ mesons, and B baryons; the lifetimes of B mesons and baryons; and the 
average B hadron decay charged multiplicity. In cc events we have considered the 
uncertainties on: the branching fraction for cc; the charmed hadron lifetimes, 

the charged multiplicity of charmed hadron decays, the production of from charmed 
hadron decays, and the fraction of charmed hadron decays containing no 7r°s. We have 
also considered the rate of production of ss in the jet fragmentation process, and the 
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Source 


Variation 


5 {xb) 


Monte Carlo statistics 




0.0011 


Tracking efficiency correction 


on/off 


0.0022 


Track impact parameter 


on/off 


0.0012 


Track polar angle 


2 mrad 


-0.0009 


Track l/P^ 


0.0017 


-0.0054 


Hadronic event selection 


standard 


0.0005 


Total Detector Systematics 




0.0061 


mass effect 


±2(T 


0.0001 


B lifetimes 


±0- 


0.0002 


B+fB^/B^Ab production 




0.0010 


B decay fraction 


±2(T 


0.0006 


B decay < rich > 


5.3±0.3 


0.0012 


D lifetimes 


±a 


0.0002 


D decay < rich > 


±a 


0.0005 


D — ^ multiplicity 


±a 


0.0013 


D no TT^ fraction 


±a 


0.0005 


g^bb 


(0.31±0.15)% 


0.0002 


g-^ cc 


(2.4±1.2)% 


0.0008 


production 


0.66±0.07 trks 


0.0009 


A production 


0.12±0.01 trks 


0.0002 


Rb 


0.2170±0.0009 


<0.0001 


Rc 


0.1733±0.0048 


<0.0001 


Model dependence 




0.0020 


Total Systematics 




0.0068 



Table 5: Source and systematic errors. 



production of secondary bb and cc from gluon splitting. The world-average values [IC, 28 



of these quantities used in our simulation, as well as the respective uncertainties, are 
listed in Table |^. Most of these variations have effect on normalisation, but very little 
on the shape or the mean value. In no case do we find a variation that changes our 
conclusion about which functions are consistent with the data. Systematic errors of the 
mean value are listed in Table ^. 

The model-dependence of the unfolding procedure is estimated by considering the 



envelope of the unfolded results illustrated in Figure [IT]. Since eight functions provide 
an acceptable probablity in ffiting to the data, in each bin of we calculated the 
average value of these eight unfolded results as well as the r.m.s. deviation. In each bin 
the average value is taken as our central value and the r.m.s. value is assigned as the 
unfolding uncertainty. 
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Other relevant systematic effects such as variation of the event selection cuts and 
the assumed B hadron mass are also found to be very small. As a cross-check, we 
vary the Momax cut (Equation (||)) in selecting the final B sample within a large range 
and repeat the analysis procedure. In each case, conclusions about the shape of the B 
energy distribution hold. In each bin, all sources of systematic uncertainty are added in 
quadrature to obtain the total systematic error. 



7. Summary and Conclusions 

We have used the excellent tracking and vertexing capabilities of SLD to reconstruct the 
energies of B hadrons in e~^e~ —>■ events over the full kinematic range by applying a 
new kinematic technique to an inclusive sample of topologically reconstructed B hadron 
decay vertices. The overall B selection efficiency of the method is 3.9%. We estimate 
the resolution on the B energy to be about 10.4% for roughly 83% of the reconstructed 
decays. The energy resolution for low energy B hadrons is significantly better than 
previous measurements. 

In order to get a good estimate of the model dependence of the unfolded distribution, 
the distribution of reconstructed scaled B hadron energy, D'^*"' [x^^^) , is compared case 
1) with predictions of ezi/ier perturbative QCD and phenomenological h quark fragmen- 
tation models in the context of the JETSET parton shower Monte Carlo, or HERWIG 
and UCLA fragmentation models, and case 2) with a set of functional forms for the B 
energy distribution. In case 1), the Lund and the Bowler models are consistent with 
the data; the model of Kartvelishvili et al. is in marginal agreement with the data. The 
models based on the perturbative QCD calculations of Braaten et al., and of Collins and 
Spiller, and the Peterson model are disfavored by the data. Although both versions of 
the HERWIG model are excluded by the data, the new version is very much improved. 
The UCLA model describes the data reasonably well. In case 2), four functional forms, 
namely the two generalised Peterson functions Fl and F2, the Peterson function, and a 
constrained 8th-order polynomial are found to be consistent with the data. 

The raw B energy distribution is then corrected for bin-to-bin migrations caused by 
the resolution of the method and for selection efficiency to derive the energy distribution 
of the weakly decaying B hadrons produced in decays. Systematic uncertainties in 
the correction have been evaluated and are found to be significantly smaller than those of 
previous direct B energy measurements. The final corrected xb distribution (^j;^"*^) 



is shown in Figure |TT[ The statistical and unfolding uncertainties are indicated sepa- 
rately. 

It is conventional to evaluate the mean of this B energy distribution, < Xb >■ 
For each of the eight functions providing a reasonable description of the data (four 
from case 1) and four from case 2)), we evaluate < > from the distribution that 
corresponds to the optimised parameters; these are listed in Table ^ and Table ^. We 
take the average of the eight values of < xb > as our central value, and define the 
model-dependent uncertainty to be the r.m.s. deviation within each bin. All detector 
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and physics modeling systematic errors are included. We obtain 

<xb> = 0.714 ± 0.005(stat.) ± 0.007(s|/st) ± 0.002(mode/), (12) 

It can be seen that < > is relatively insensitive to the variety of allowed forms of the 
shape of the fragmentation function D{xb)- 
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Figure 1: Distribution of the reconstructed PfCorrected vertex mass in the 1996-97 
data (points). Also shown is the prediction of the Monte Carlo simulation, for which 
the flavor composition is indicated. 
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Figure 2: The relative deviation of the maximum missing mass from the true missing 
mass for Monte Carlo simulated B hadron decays, which is divided into three categories: 
and B"^ (open), B^ (cross-hatched), and A;, (dark hatched). 
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Figure 3: Distribution of the reconstructed Mq^^^, for the selected vertices in the 
1996-97 data (points). Also shown is the prediction of the Monte Carlo simulation. 
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Figure 4: Distribution of the reconstructed Mq^^^ for the final selected B sample (see 
text). Also shown is the prediction of the Monte Carlo simulation. 
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Figure 5: The Monte Carlo simulated efficiency for selecting B hadron decay vertices 
as a function of the true scaled B hadron energy, x*™*^ = E*j^^^ / Ebeam- The nearly 
energy-independent efficiency (except at very low B energy) improves the sensitivity of 
the measured x^'^ distribution to the true underlying B energy distribution. 
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Figure 6: The fitted core and tail widths of the B energy resolution as a function of the 
true scaled B hadron energy. The ratio of the amplitude of the inner Gaussian (core) to 
that of the outer Gaussian (tail) is 83:17. The dependence of the core resolution on the 
true B energy is small. The very good resolution for low energy B hadrons improves the 
sensitivity of the measured a;^'^ distribution to the true underlying B energy distribution. 
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Figure 7: Distribution of the reconstructed scaled B hadron energy for 1996-97 data 
(points) and the default Monte Carlo simulation (histogram). The solid histogram shows 
the non-66 background. 
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Figure 8: Each figure shows the background-subtracted distribution of reconstructed B 
hadron energy for the data (points) and for the Monte Carlo (histogram) based on the 
respective optimised input fragmentation function within the JETSET parton shower 
simulation, as well as based on the HERWIG {dd = and eld = 1) and the UCLA 
fragmentation models. The ^-^id the number of degrees of freedom are indicated. 
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Figure 9: Each figure shows the background-subtracted distribution of reconstructed B 
hadron energy for the data (points) and for the weighted simulation (histograms) based 
on the respective optimised input functional form for the true B energy distribution. 
The the number of degrees of freedom are indicated. 
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Figure 10: The efficiency-resolution corrected distributions of scaled weakly- decaying 
B hadron energies for Case 1) fragmentation models of the Lund, the Bowler and the 
Karvelishvili within the JETSET parton shower Monte Carlo as well as for the UCLA 
fragmentation model; and for Case 2) four functional forms: Fl, F2, Peterson, and the 
constrained 8th-order polynomial. 
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Figure 11: Distribution of the final corrected scaled B hadron energies. The central 
value is the bin-by-bin average of the eight consistent B energy distributions. In each bin 
the statistical error is indicated by the inner error bar, the sum in quadrature of statistical 
and unfolding errors from model dependence by the outer error bar. Systematic errors are 
small compared with the statistical and model dependence errors and are not included 
here. Note that the first two bins are below the kinematic limit for xb (no point shown). 
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